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ABSTRACT 

We describe a practical measurement of the curvature of the Universe which, un¬ 
like current constraints, relies purely on the properties of the Robertson-Walker metric 
rather than any assumed model for the dynamics and content of the Universe. The ob¬ 
servable quantity is the cross-correlation between foreground mass and gravitational 
shear of background galaxies, which depends upon the angular diameter distances 
dA{zi), dA{zs), and dA{zs,Zi) on the degenerate triangle formed by observer, source, 
and lens. In a flat Universe, dA{z£, Zg) = — dA{z£), but in curved Universes an 

additional term oc 12^ appears and alters the leasing observables even if dA{z) is fixed. 
We describe a method whereby weak leasing data may be used to solve simultaneously 
for Oa and the curvature. This method is completely insensitive to the equation of state 
of the contents of the Universe, or amendments to General Relativity that alter the 
gravitational deflection of light or the growth of structure. The curvature estimate is 
also independent of biases in the photometric redshift scale. This measurement is shown 
to be subject to a degeneracy among dA-, flfc and the galaxy bias factors that may be 
broken by using the same imaging data to measure the angular scale of baryon acoustic 
oscillations. Simplified estimates of the accuracy attainable by this method indicate that 
ambitious weak-lensing -|- baryon-oscillation surveys would measure to an accuracy 
~ 0-044ky'^^(crin^/0.04)^/2^ where ainz is the photometric redshift error. The Fisher- 
matrix formalism developed here is also useful for predicting bounds on curvature and 
other characteristics of parametric dark-energy models. We forecast some representative 
error levels and compare to other analyses of the weak leasing cross-correlation method. 
We find both curvature and parametric constraints to be surprisingly insensitive to the 
systematic shear calibration errors. 

Subject headings: gravitational leasing; cosmological parameters; relativity 
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1. Metric Measurements of Curvature 


The most robust prediction of inflation theories is that the radius of curvature of the present 
Universe should be very large, i.e. = 0 to high accuracy. The WMAP+supernovae+i^o best 
fit of = —0.02 ± 0.02 (Spergel et al. 2003) is therefore taken as a vindication of inflation, and a 
great majority of current work, such as investigations of the properties of dark energy, take = 0 
as a priori truth. It is important to realize, however, that WMAP does not measure curvature 
in any direct geometric way. Constraints on curvature derive primarily from a measurement of 
the angular diameter distance to recombination, DA{zrec), which depends upon curvature but also 
upon Qm and models for dark energy or other constituents of the Universe. Hence the curvature 
results are dependent upon the ACDM (or other) model for dark energy assumed in the analysis. 
Our ignorance of the dark energy phenomenon limits our ability to test for flatness or to look for 
small finite that may be predicted in variants of inflation (Uzan, Kirchner, & Ellis 2003) or in 
Universes with non-trivial topology (Luminet et al. 2003). 

Our ignorance of the true curvature will, conversely, foil attempts to characterize dark energy 
properties with Type la supernovae (or other standard candles), which measure Dl{z) at lower 
redshifts. If the Tolman surface-brightness relation holds, then we have Da = Dl{1 + The 

proper-angular-diameter distance in a Universe with Robertson-Walker metric is 
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Here r is the comoving radial distance, Rq is the radius of curvature of the Universe, and H{z) = 
a/a. These equations follow purely from the RW metric. If standard-candle and CMB observations 
were to give us perfect knowledge of Da{z), then for any Rq and k there is a solution 

^ = + (4) 


that will exactly reproduce the data. If the Friedmann equations hold, then 
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Here fxiz) describes the evolution of px, an additional dark-energy component. Hence even with 
perfect knowledge of Da{z), there is always some dark energy behavior fx which reproduces that 
data for any choice of curvature. The degeneracy between curvature and dynamics is well known, 
e.g. Weinberg (1970) demonstrates that even complete knowledge of Da{z) cannot constrain the 
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curvature in the absence of a dynamical model for the expansion. Constraints on curvature from 
Da data arise solely because of our preferences for forms of fx that arise from certain equations 
of state. Allowing alterations to the Friedmann equations adds additional degeneracy. 


We can break the degeneracy between curvature and fx without making dynamical assump¬ 
tions if we apply the metric to some line segment that does not originate at z = 0, essentially 
forming a cosmological-scale triangle. For example if a gravitational lens at Z£ bends a light ray by 
angle a, then the apparent deflection of a source at Zg observed from z = 0 obeys 


69 = a 


DA{zi,Zs) 

Da{Zs) 


(7) 


where Da{z^,Zs) is the proper angular-diameter distance at Zg as viewed from z^. We introduce a 
dimensionless comoving angular diameter distance 
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Now xo is the radius of curvature in Hubble lengths. Using the difference rule for S{x) we obtain 
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So while the d{z) function itself has a curvature-dark energy degeneracy, the lensing strength 
depends purely on once d{z) is known, e.g. from supernova measurements. This dependence 
is well known, e.g. Equation (13.70) in Peebles (1993). Figure 1 illustrates how the curvature of the 
Universe is indeterminate when distances from z = 0 are the only observables, whereas a measure 
of dig determines curvature. Linder (1988) remarks that the presence of this differential distance 
in gravitational lensing equations could provide useful cosmological tests, but to our knowledge no 
practical implementation of this—or any other—metric curvature test has been proposed. 
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Fig. 1.— Illustration of curvature determination: consider an observer on a section of the closed Universe, 
represented by the bold line. The conformal distances xl and xs to two sources are not directly observable; 
only the angular-diameter distances are, cIl and ds represented here by the chord lengths. The radius of 
curvature xo is indeterminate, but a measure of the distance (dashed line) allows determination of the 
curvature. 

We propose that be determined by measuring the gravitational leasing shear 7 of back¬ 
ground galaxies— i.e. the gradient of 66 in Equation (7)—and measuring the amplitude of the 
contribution from the O(nfc) term in Equation (17). A foolishly optimistic estimate of the ac¬ 
curacy available on can be obtained by assuming that the deflection angles a are known a 
priori at all redshifts and positions on the sky, and that the angular diameter distance d{z) has 
been measured as well. Since dg and dg are in the range 0.5-1.5 at relevant distances for weak 
leasing in the ACDM fiducial model, the lensing shear scales as ~ (1 — 14^/2), so we would find 
that 6Qk ~ 257 / 7 . The typical shear amplitude 7 over a line of sight is ~ 0.02 at cosmological 
distances. The uncertainty in shear is « 0.3/^/Ng, where Ng is the number of galaxies with well- 
measured shapes. A space-based survey can measure 100 galaxies per square arcminute, which leads 
to 6Qk ~ 0-0002/-v/7g]^y, far better than the current model-dependent constraints. This estimate 
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is unrealistic, however, because neither the distance functions nor the deflections will be known. 
Below we will investigate a simultaneous solution for the angular-diameter distances, deflection 
strengths, and flfc from cross-correlation of lensing shear with foreground galaxy distributions, and 
estimate the resultant degradation in the curvature constraint. 

If we execute the WL curvature measurement by dividing our galaxy sample into redshift bins, 
then we note that we do not actually need to know the redshift of each bin; it is only the distance 
d of each bin that enters into the calculation of the shear. Hence biases in photometric redshifts 
do not affect the curvature measurement. It is essential, however, to insure that the galaxies in a 
bin are at a common distance. 

In the remainder of this section we examine metric tests for curvature which might be possible 
with other cosmological observables besides weak lensing. In §2 we develop the methodology for 
using weak lensing cross-correlations to solve simultaneously for the curvature, distance relation¬ 
ships, and deflector properties. The implementation is significantly more involved and subtle than 
the simple idea that curvature is manifested in d^g- §3 applies the formalism to forecast curvature 
constraints from feasible surveys, and extends the formalism to include treatment of the likely 
dominant systematic error. §4 investigates the constraints on parametric models of dark energy 
{wo,Wa) that can be derived with the curvature left free to vary, and compares the present results 
to some previous work on weak lensing cross-correlations. §5 summarizes and concludes. 


1.1. Metric Curvature from Other Cosmological Measurements 

Similarly metric determinations of are in principle derivable from other cosmological mea¬ 
surements. The most promising is the detection of the recombination acoustic horizon scale, which is 
known in comoving physical units from the CMB, in the power spectrum of galaxies, recently demon¬ 
strated by Eisenstein et al. (2005). The transverse baryon acoustic oscillation (BAO) test measures 
the angular scale subtended by this standard ruler, and hence measures (1 -|- z)Da{z) oc h~^d{z). 
The transverse BAO method does not by itself break the curvature-dark energy degeneracy, but 
may serve as a source of high-reliability constraints on d{z) to combine with weak lensing infor¬ 
mation. Again for a given galaxy subsample, the redshift itself is unimportant, just the observable 
d. 


With sufficient redshift resolution, the acoustic scale may be identified along the line of sight, 
yielding knowledge of 

^ = -^^d'{z){l - nkdV2) + 0(H|). (18) 

Hence the comparison of line-of-sight to transverse BAO scales can yield directly, but requires 
one to either integrate the line-of-sight H(z) data or differentiate the transverse d{z) data. Doing 
either operation in the presence of noise adds substantial difficulties to making a high-accuracy 
comparison, unless one assumes a parametric form for d{z). 
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Counts of galaxy clusters are sensitive to dV/clQ dz oc h~^d‘^d'{z){l — Qkd‘^/2). Volume-element 
data alone can not measure 17^, but could potentially do so in combination with high-precision 
measures of d{z) and its derivative. Cluster counts are also, however, dependent upon the growth 
of structure, so do not offer the model-independent cosmographic information of weak leasing or 
baryon oscillations. 

Stellar evolution has been proposed as a standard chronometer (Jimenez et al. 2003), which 
could determine H{z) and hence the curvature if combined with measures of d. To contribute to 
the metric measurement of curvature, accuracies of a few percent on H{z) are required, and further 
investigation into the complexities of stellar and galactic evolution are needed in order to determine 
whether robust constraints at this accuracy are possible. 

Strong lensing may also be used to discern the behavior of d{zt, Zs)/d{zs)., if one is presented 
with a lensing system with multiply-imaged sources at a variety of (known) background redshifts. 
Such use of multiple arcs around clusters has been investigated in the context of constraining dark 
energy (Link & Pierce 1998; Gautret, Fort, & Mellier 2000; Golse, Kneib, & Soucail 2002; Sereno 
2002 ; Soucail et al. 2004), and would be equally applicable to determination of curvature, but also 
equally susceptible to the small uncertainties in the mass profiles of clusters. 


2. Cosmographic Methodology 

The metric curvature determination is clearly an extension of the cross-correlation cosmog¬ 
raphy (GGG) technique proposed by Jain & Taylor (2003): one identifies rich foreground clusters 
and measures the dependence of the induced background shear upon of the source galaxies. 
The mass of the cluster(s) may be unknown, but the ratio of any two background shears is a 
purely cosmographic function. Song & Knox (2004) combine this approach with other dark energy 
probes. Bernstein & Jain (2004) [BJ04] generalize the method to a cross-correlation between fore¬ 
ground estimated-mass distributions and background shear patterns. In this case the cluster-mass 
uncertainty is replaced by an unknown bias factor for each foreground mass shell, over which we 
marginalize. Zhang, Hui, & Stebbins (2003) [ZHS] give a more elegant analysis in which the galaxy- 
shear cross-correlations are considered simultaneously with the shear-shear correlations in a single 
Fisher matrix. Hu & Jain (2004) also include the galaxy-galaxy correlations in the same Fisher 
matrix, and make use of a parametric model of bias based on the halo model of galaxy distributions. 
All of these analyses differ substantially in their analytic approaches, underlying assumptions, and 
estimated constraints. A full comparison is beyond the scope of this paper, but we will make 
some brief comments below. Here we attempt to recast the problem in a manner suited to the 
metric curvature constraint, but also amenable to studying the parametric dark-energy constraints 
investigated by these authors. 

The shear induced in direction 6 on sources at Zg by mass in the foreground can be decomposed 
into E- and H-mode components. In the weak-lensing limit, the former is straightforwardly related 
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to the lensing convergence k and the latter should vanish. Hence we will consider the shape 
information to come in the form of the convergence field k, which depends upon direction and 
source redshift as 


K{e,zs) = 


j d(,{l + Z()Ho ds-di ^ ^ 

dze -f 77 -i- rn{6, ze )—--(1 - nkdsde/2), 


H{ze) 


do 


(19) 


where m is the overdensity 6p/p of the total mass. 

We also idealize the sky as consisting of K discrete pencil-beams of solid angle dH, with the 
true mass density m{6i) being statistically independent of the density in any other beam. We can 
choose dH to correspond to a correlation length of the true continuum shear field. We will assume 
dH = 0.25 arcmin^; the results depend only weakly on this choice. We could equivalently decompose 
the shear and mass distributions into spherical harmonics, in which case our choice of 5H becomes 
equivalent to a maximum multipole i ^ 2 x 10^. 

The weak lensing equation (19) can be written 


K{s,6i) = J d£G{s,£)m{£,6i)k{£), 

( 20 ) 

k{£) = ^£{1 + z£){i - nkeV2) 

( 21 ) 

G{s,£) = {l-nks£/2)eis-£), 

( 22 ) 


where 0 is the step function, and we have reparameterized the source- and lens-plane distances by 
s = d{zs), i = d{z£). This equation assumes the validity of General Relativity in two respects: first, 
that light follows the geodesic equation for the metric, and second, that we have the usual Poisson 
equation. The second assumption can be dropped in our quest for a purely metric test of curvature. 
In our formalism a change of the Poisson equation would be accommodated by generalizing k£m{i, 0) 
to mean the convergence caused by the mass at distance i, even if this means that m is no longer 
the unaltered mass distribution. 

We discretize the analysis by considering the galaxies to be located on a series of J shells, 
in order to facilitate the construction of Fisher matrices below. The finite thickness of the shells 
complicates the analysis (ZHS) and in fact it is possible to draw erroneous conclusions from a 
discrete analysis, as discussed by Stebbins (2005). Below we will attempt to return the discrete 
formulation to the continuum limit and discuss these issues. 

The measured shear at each source plane becomes, with a noise term now added, 

Gs£m£{0i)ki + {6K)s{0i), (23) 
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\/n75rT’ 

ke = ^^de{l +ze)Ax = ^^de{l +ze){l-0.kdj/2)Ade. (26) 


Now mi{6) is the mass overdensity averaged through the shell of width Ax; k(, is a factor that 
converts this overdensity into the convergence measured on a lens plane at infinity. The shear 
measurement noise is determined by the areal density Ug of shape-measurable galaxies on shell s, 
and cr^ 0.3 (Bernstein 2005). 

The true mass fluctuation m(,{0) is not directly observable. But we assume that we deduce 
from the galaxy distribution some estimate of the mass overdensity ge{0) that is correlated with 
the true distribution. This proxy field need not be the galaxy density itself, but rather the result 
of any algorithm applied to the observed galaxy field, e.g. involving the assignment of halos to 
galaxies, groups, and clusters. The correlation coefficient of the true and proxy mass fields is 


imgj) 


— dijfi, 


(27) 


where the angle brackets denote averaging over direction 6 . The r* are not known but are assumed 
positive. We assume that the shells are thick enough to be uncorrelated. We also define for each 
shell an unknown bias via 


{ruigi) = Bi{gf). 


(28) 


Note that this is the bias of the true mass with respect to the estimated mass. The more conventional 
bias b of the estimated (galaxy) mass to the true mass is, for Gaussian-distributed g and m, related 
to this via 



(29) 


The shear noise 6k is assumed to have correlations with neither the true nor the proxy mass 
distribution. 


With these definitions we could extract curvature information from our observed shear data 
as follows: we cross-correlate each convergence held Kg with each foreground proxy held gi to give 
an observable quantity Xg£. From Equation (23) the expectation values are 

i^si) = Ggiikii^gimji) = GgiBiki{gj). (30) 

l' 

There are J(J — l)/2 non-zero observable cross-correlations, which are statistically independent for 
any signihcant sky coverage. These Xgi are determined by the J factors {B£ki{g‘j)), plus the matrix 
G that is dehned by J distances di and the single hlfc. Hence there are 2J -|- 1 free parameters, so 
we might expect an unambiguous ht to the data for J > 5. After marginalization over the nuisance 
parameters Bi and di (which carry much information about dark energy and galaxy formation!) 
one obtains constraints on 
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Unfortunately the simultaneous solution for 17^, B, and d from the system of equations (30) 
has three degeneracies in the neighborhood of = 0. The transformations 

di^di{l + ao) (31) 

di ^ di/{l - aidi), {Biki) ^ - aidi) (32) 

di ^ di/{l - a2d?i), {Biki) ^ {Biki){l - a2d'f), 17^ ^ flfc + 2a2 (33) 


each leave the observable shear unchanged to first order in 17^ and the free parameters a*. The first 
degeneracy is a simple scaling of the di, which is unsurprising as the leasing observables depend 
only upon distance ratios. The second degeneracy is similarly benign as it leaves the solution for 
f7fc unchanged. But the third degeneracy leaves 17^ indeterminate. The degeneracy can be broken 
if we assume some functional form for d{z), but our goal is a model-independent constraint. The 
same imaging survey that is used to generate the shear measurement and the foreground galaxy 
maps can be used to measure the transverse baryon acoustic horizon scale, if the photometric 
redshift accuracy ai^z is sufficiently good. This will produce a statistically independent measure of 
di at each shell, which we may use to break the curvature degeneracy that remains in the lensing 
cross-correlations. 

There is a more subtle problem with a solution for cosmology via Equation (30). As noted by 
Stebbins (2005), the solution offers no discriminatory power on cosmology if we allow the Bi to be 
completely free, because the matrix Gge must be invertible as we approach the continuum limit. 
Looking explicitly at the continuum limit in Equation (20), Stebbins shows that there must always 
be some function s) such that 


m{i)k{£) = dsG ^{i,s)K,{s). 


In fact we can construct this function explicitly to first order in 17^; 

G-\£, s) = (^1 + ^£^'^ s5"{s -£) + VLk£s5'{s - £) 

d{sK 


f7i. d?‘{sK) 


m{£)k{£) = ( 1 + 


ds"^ 




ds 


(34) 


(35) 

(36) 


Here 6 is the Dirac function. Surprisingly this inversion of the lensing cross-correlation is local. 
Given this solution, we simply set B{£) = m{£)lg{£) and get an exact solution on this line of 
sight, regardless of our choice of 17^ and d{z) functions. There is hence no possibility of inference 
of cosmological parameters along this line of sight as we reach the continuum limit. Similarly a 
solution of the cross-correlation Equation (30) is possible for any choice of cosmology, once we reach 
a large number of lens and source planes, as G is always invertible to give a valid B. 


An exit from this conundrum is the realization that giving complete freedom to the bias 
function B{£) is tantamount to discarding our initial presumption that the proxy field is correlated 
with the true mass distribution. Imagine slicing the Universe so finely as to assign each proton 
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its own bias value with respect to the dark matter. With this freedom we could clearly produce 
any macroscopic mass distribution we wish, including one that has no correlation at all with the 
galaxy-scale baryon distribution. We must therefore incorporate into the analysis some criteria for 
the coherence of the bias parameter in order to reflect our underlying assumption that mass has 
some correlation with the proxy field. 

We note, for example, that the continuum-limit inversion formula Equation (35) will produce a 
divergent solution for m{i) in the presence of any finite amount of shot noise on the lensing field k{s). 
The cross-correlation between estimated m{£) and g{£) can be estimated from the multiple lines of 
sight in the survey, and will be driven to zero in this continuum limit. If the likelihood function 
for our joint analysis of the k and g fields requires them to be correlated Gaussian variables, then 
the rapidly varying components of the inversion solution for m will be suppressed in a maximum- 
likelihood solution, and it becomes possible to discriminate cosmologies. 

We do not claim yet to have a rigorous model-independent method to approach the continuum 
limit. Below we produce a likelihood function for the shear and proxy fields that incorporates their 
presumed correlation on scales of a shell thickness. Our practical approach will be to see how the 
uncertainties in our cosmological parameters scale as the number of shells increases; in fact we see 
no significant changes in any of the results below as we decrease the shell size as far as Az ~ 0.02. 
Another practical approach that we will take is to include a regularization condition on B{i) to 
lower the probability of solutions that vary rapidly on short redshift (and hence time) scales, in 
accordance with our physical intuition; this too has no significant impact upon our forecasts. 

A potentially more enlightening approach to the cross-correlation problem may be to param¬ 
eterize the functions d{z) and B[z) by expansions in orthonormal function sets, such as Fourier 
modes or Legendre polynomials, rather than assume a stepwise variation and/or discrete galaxy 
distributions. 


2.1. Likelihood Function 

The accuracy of constraints on can be estimated by the Fisher matrix methodology. We 
produce a Fisher matrix for a single line of sight, and then multiply by K since we have assumed 
all lines of sight to be independent. On each line of sight there are three J-dimensional vectors: 
the measured convergence k] the true mass distribution m; and the proxy mass estimate g. We 
need a likelihood function L{k, in,g), which we then must marginalize over the unobservable m. 
This probability can be expressed as 


L{K,ui,g) = L(K:|m)L(m,g). 


(37) 
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Since the shear measurements arise from the sum of many individual galaxy shape measurements, 
the shear measurement noise is well expressed as a Gaussian distribution, so the first term is 

L{K\m) = (27r)-‘^/2|N|"^/2exp|-i [(k - Gkm)'^N“i(K - Gkm)] | . (38) 

where N is the convergence noise matrix diag(((5K^)), and k = (X\a.g{k(). 

For analytical convenience we take the joint probability L(m, g) of the true and estimated 
masses to be a bivariate Gaussian at each slice. This is undoubtedly a poor approximation, e.g. 
most of the lensing cross-correlation power is on angular scales where the mass distribution is 
highly non-linear and non-Gaussian. We defer a more realistic treatment to later work. Given our 
definition of the bias Bi and correlation coefficient r*, the probability can be written 


L(m, g) = (27r)-'^|Cur^/"|Cg|-i/2 exp |-i [(m - C^-\Tn - Bg) + g^Cg'^g] |(39) 

Cg = dmg{{gj)) (40) 

B = diag(Bi) (41) 

Cu = diag = diag [{mj){l - rf)] . (42) 


Cg is the correlation matrix for the mass estimator, and Cu is the correlation matrix for the part 
of the mass fluctuations that are uncorrelated with the estimator g. We assume the shells are thick 
enough for both matrices to be diagonal. 


Multiplying Equations (38) and (39) , then integrating over m gives the Gaussian distribution 

L(K,g) = (27r)“‘^|KCg|~^/^exp |-i [g'^Cg"^g-F (k - GkBg)'^K~^(fv - GkBg)] I , (43) 
K = N + GkCuk^G^. (44) 


K is the covariance matrix for k if the mass estimators g are held fixed, i.e. only the components 
of the mass distribution that are uncorrelated with g are considered stochastic. 


From this multivariate Gaussian distribution we may derive a Fisher matrix using the formulae 
for zero-mean distribution given, e.g. , by Tegmark, Taylor, & Heavens (1997). We also multiply 
this by the number of independent lines of sight in the survey to get 


F- • — 


47r/sky 

sn 


Tr 


Cg-'Cg^,C, 


- 1 / 


■'s,j 


+ K"^K 


+ 2K~^(GkB),iCg(GkB) 


(45) 


with the commas in the subscripts denoting differentiation. The Fisher information nicely separates 
into three parts; the first is the information that can be gleaned from the variances of the mass 
estimator g, i.e. the galaxy power spectrum. The third term is Fisher information that would 
arise from adjusting parameters to minimize the the fit of k to the estimated mass g in 

Equation (23), were the values of g taken as fixed and the matrix K taken as a known covariance 
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for the K values. The third term is information gleaned from the actual covariance of the k residuals 
to this fit, and looks just like the Fisher matrix for pure shear power-spectrum tomography, except 
that the relevant mass power spectrum in this term is Cu, the power that is uncorrelated with the 
galaxies, not the (larger) total power of m. Hence the inclusion of cross-correlation will lower the 
sample variance in lensing power-spectrum tomography. 

In this paper we will presume that none of the (co-)variances of the true or estimated mass 
distributions are well predicted by theory. Hence we will marginalize over Cg,B, and Cu at each 
redshift. There will be no cosmological information gleaned from analysis of Cg since we consider 
its elements to be free parameters of the model, hence we can drop the first term of Equation (45). 

Since we are marginalizing over all the mass variances, we can also absorb the overdensity-to- 
convergence conversion factor ki into the definition of m£ and g£, and remove it from Equations (44) 
and (45). We just must remember to express our fiducial values of Cu and Cg as convergence 
variances rather than as over density variances. 

In a future paper we will consider cases where there are reliable theoretical models for the 
mass fluctuations. Hu & Jain (2004) consider the case where not only the mass power spectrum 
but also the bias and correlation of the galaxy sample are described by parametric models. 

To recap, then, we will use the Fisher matrix 

Fij = ^^^Tr [K-1 k,K-1k,+ 2K-i(GB),iCg(GB);^] , (46) 

K = N + GCuG^. (47) 

with the parameters being; 

• Hfc, which affects only G. The fiducial value is zero. 

• The di, which also only affect G. The fiducial values at each redshift are taken from the 
ACDM model with = 0.3. 

• The diagonal elements of Cu, over which we will marginalize to allow for our ignorance of 
the cross-correlations r£ between the proxy and true mass distributions. We take a fiducial 
r£ = 0.7 but investigate the effect of different fiducial values. 

• The diagonal elements Bi of B, the bias factors, over which we will marginalize. We take the 
fiducial Bi = 1 without loss of generality. 

• The diagonal elements of Cg are also free parameters, but striking the first term of (45) is 
equivalent to marginalizing over these parameters. 


With these parameters, all of the derivatives required in Equation (46) are simple, sparse matrices. 
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2.2. Additional Information 


The CCC Fisher matrix will be singular due to the degeneracies in the parameters noted above. 
To this we must add a Fisher matrix for the di determined from the baryon acoustic oscillations; 
this matrix will be diagonal as the galaxy power spectrum measurements are independent from 
shell to shell if the shells are much thicker than the acoustic horizon scale. 


For the baryon acoustic scale measurement, we use an abstraction of the full analysis provided 
by Seo & Eisenstein (2003). We presume that there is a fractional uncertainty (Tin(i+ 2 ) in 1 + z 
on each photometric redshift, which defines the depth of a “slab” of modes in /c-space that are 
measurable. We also presume that the power spectrum determination for all modes within the slab 
will be limited by sample variance, which requires that the comoving volume density of galaxies with 
good photo-z’s must be dN/dV > Mpc“^ = to make the shot noise negligible at 

the highest k ~ 0.2/iMpc~^ where baryon oscillations are usefully detected. When these conditions 
are met, the fractional error in the mean d for a given shell will scale as \J(y\n{i+z)/^■ The prefactor 
depends upon the intensity of the wiggles in P{k), which diminishes somewhat at lower redshift as 
nonlinear effects erase the acoustic features. We approximate this by assigning an uncertainty 

for those shells which, by Equation (50), have dN/dV above the sample-variance-limited threshold 
noted above. If the galaxy density in the shell is below this level, we assume no constraint on d 
from BAO. Our simplification fits the Seo & Eisenstein (2003) results well and comes within 30% 
of the forecasts by Glazebrook & Blake (2005). 


We will also investigate the impact of a regularization constraint on the B function. Writing 
the bias as a function of the lens-plane angular diameter distance i, we add to the likelihood a 
function 


1 


The parameter R specifies a scale of the RMS bias slope that will be suppressed. The discretized 
version of this likelihood is a quadratic form in the Bi and easily incorporated into the Fisher matrix. 
We project the three degeneracies in Equations (31-33) out of the B vector before applying the 
regularization constraint, so that the regularization does not artificially break these degeneracies. 


3. Curvature Forecasts 

We calculate the Fisher matrix about our fiducial cosmology using the following estimates for 
the relevant quantities. 
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The K variances in N are taken from Equation (25) with the source galaxies distributed as 

^ oc 2 :^ exp [-{z/lAlzof-^] , (50) 

with the median redshift zq and the integrated density n chosen to crudely mimic the expectations of 
future surveys. Our default is to take zq = 1.5, n = 100 arcmin“^, the performance one might expect 
from a space-based survey like the proposed SNAP wide survey (Massey et al. 2004; Bernstein 2005). 

We set the variance (m|) to the convergence variance d{K^)/dz integrated over each shell, using 
the nonlinear dark-matter power spectrum for a ACDM model, restricted to A: < 27r/(50 kpc), as 
detailed in BJ04. Given a fiducial correlation coefficient r, we set ((/|) = We will set all 

the fiducial r£ = 0.7, although the marginalization allows them to have independent solutions. 

We also truncate the source and lens distributions at Zmax = 3 unless otherwise noted. The 
number of redshift bins will range from 30 to ~ 200 to search for effects of the artificial discretization 
of the problem. 


3.1. Constraints on 

The Fisher matrix for the fiducial cosmology yields, first, an estimate of the uncertainty on 
flfc of cjfc = 1.6 X in the unrealistic case where all the biases, distances, and correlation 

coefficients are known a priori. This is in good agreement with our crude guess earlier. 

When we marginalize over all parameters except 14^) we find that if we project out the degen¬ 
eracy in Equation (33), then the weak leasing CCC method is remarkably efficient at extracting 
flfc) giving an uncertainty cr^ ~ 2 x nearly as good as before any marginalization. Of 

course this is akin to stating that the Titanic was watertight, aside from the big gash from the 
iceberg. Nonetheless this feature of the CCC method is illuminating: it tells us that the role of 
other distance measures, i.e. supernovae and transverse BAO, will be purely to bound the size 
of the quadratic degeneracy in the CCC distances that couples to 14^, because the CCC method 
reduces any other form of uncertainty in d{z) and 0^ to far below the level achievable from such 
experiments. 

If, for example, we remove the three singular values from the WL-only Fisher matrix, and then 
estimate the uncertainties in d{z), we find ~ 1.4 x over the range 0.6 A z < 2.5, when 

the distances are averaged over Az = 0.1 bins. This is 1-2 orders of magnitude better than the 
most ambitious targets for supernova measurements, and ~ lOx better than a spectroscopic BAO 
survey of equivalent sky coverage. 

The final constraint on 14^, therefore, comes down to the simple question of how well the 
non-WL data can constrain a variation of the form 


d —> (i(l -\- cxq -\- cx\d -\- r4/;,(i^/2). 


(51) 
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with the WL CCC data having locked down any other form of variation. Our simplified analysis 
of the transverse BAO survey available from the canonical WL survey data yields 



In our canonical survey, both WL information and sample-variance-limited transverse BAO 
data are available over the range 0 < 2: < 3. In Figure 2 we examine the dependence of the 0,^ 
constraint on the redshift range over which these data are jointly available. Reaching a 1% model- 
independent constraint on 0^ would require a massive spectroscopic BAO survey, and/or obtaining 
CCC-I-BAO data to 2 ~ 8. This is not a completely ridiculous prospect: the CMB anisotropy is, 
after all, a transverse BAO measurement at 2; = 1100. Leasing of the CMB is also observable, but 
even more leasing information may be available using the 21-cm emission from neutral hydrogen 
in the early stages of reionization (Pen 2004). 

Figure 2 also plots the 0^ limits that result from combining the transverse BAO distance 
measurements with a supernova survey that yields independent 1% measurements of d{z) at A2; = 
0.1 steps in the range 0.4 < 1.7. We find that the SN data can cut cTk in half for a transverse BAO 
survey with /g^y = 0.1, but the SN data are of only slight help (20%) to a full-sky photo-2; BAO 
survey. 

We note that the Fisher uncertainties on 0^ are quite independent of the number of redshift 
shells or of any regularization constraint on the bias factor. They are also essentially independent 
of the assumed proxy correlation coefficient r or the WL shape noise a-y, because the nominal WL 
survey is already reducing the errors to insignificance in any mode for which it has power at all. 
The galaxy distribution parameters zq and n are important only in that they have been used to 
determine the redshift range over which sample-variance-limited BAO data are available. If we scale 
back to zq = 1.0, n = 40arcmin“^, as might be expected for a deep ground-based survey (Bernstein 
2005), then the uncertainty in 0^ triples to 0.12/^^,^ for the nominal a\nz = 0-04 because the galaxy 
density drops below the BAO sample-variance limit at 2; < 3. In a real survey, the distribution 
of galaxies with accurately measured shapes will not, however, equal the distribution of those 
with accurately measured photometric redshifts, so a more detailed analysis is required in order to 
compare surveys. It is clear, however, that sky coverage, photo-2; accuracy, and redshift range of 
joint CCC-BAO data are the key characteristics for this geometric constraint of 0^. 


3.2. Effects of Systematic Errors 

The WL CCC method depends upon measuring shear to very high accuracy. Indeed the 
cosmology dependence of shear is so subtle that Mandelbaum et al. (2005) reverse the technique, 
and use the redshift dependence of galaxy-shear correlations as a test for systematic errors in shear 
measurement. Any such shear systematics could substantially degrade the curvature constraints 
derived above. But we show here that the CCC data will be rich enough to solve simultaneously 
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Fig. 2.— Uncertainties obtained on are plotted vs the maximum redshift for which weak lensing cross¬ 
correlation and baryon acoustic oscillation data are jointly obtained. These are purely metric constraints 
on Ufe, and it is presumed that the WL data constrain all forms of variation of except the principle 
degeneracy described in the text. The upper two curves are for 10% sky coverage, the lower two for full-sky. 
In each pair the solid line is for WL-I-BAO constraints, and the dashed line shows the effect of adding 1% 
measures of supernova distances in each /S.z = 0.1 bin for 0.4 < z < 1.7. 


for cosmology and calibration systematics. 

The WL cross-correlation technique is insensitive to spurious shear signals caused by PSF 
ellipticities, because they will not correlate with the proxy mass helds. Calibration errors on the 
shear will, however, be important. We examine a model, similar to that of Ishak et al. (2004), in 
which the shear on source shell s is mismeasured by a factor {1 + fg). This is simply incorporated 
into our governing equations by setting 

Ggi = {l + fs)^^^{l-nkdsdi/2) (53) 

CLs 

for dg > di- The Fisher matrix now has the parameters /* as well. We assume a Gaussian prior 
distribution on the /* for which they are all independent, with Var/j = a'j. 
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We find the constraints on 11^ from the canonical survey are essentially unchanged if cjj < 
0.01/gky' , a level that is probably achievable from ground or space observations. For uj = 1, i.e. 
complete ignorance of the calibration, ak rises to 0.067 for the full-sky survey. The degradation 
in the curvature constraint is mild because even a “self-calibrated” CCC survey constrains the 
non-degenerate modes fairly well, and the nature of the degeneracies is only slightly expanded by 
the additional free parameters. A strictly z^-dependent calibration error does not appear to be a 
problem for the curvature measurement. 


4. Parametric Modelling 

The Fisher matrix constructed over parameters di and marginalized over the nuisance 
parameters, fully describes the cosmological information content of a CCC survey. Distance in¬ 
formation from transverse BAO and supernovae can be summed into the same matrix. If we now 
consider each di to correspond to a known redshift Zi (which was unimportant to the metric curva¬ 
ture determination), this Fisher matrix describes our knowledge of the expansion history, and can 
be used to forecast constraints on parametric models of the expansion of the Universe. To project 
the d~Qk Fisher matrix onto some parameters pi, we simply need to know dDA{z)/dpi. Calibration 
errors on the shear can be included as described above; biases Azp{z) in the photometric redshift 
scale are also easily included as they give rise to perturbations in the measured distances equal to 
[d{di)/dz] Azp{z). 

The distance-based analysis is also useful in that it has revealed the nature of the constraints 
from the various methods. The WL cross-correlations are completely blind to changes in log d that 
are quadratic (linear) in d when is free (fixed), but offer very strong limits on any other forms 
of variation in d{z). When WL is used in combination with transverse BAO and/or supernova 
measurements, the latter will serve primarily to constrain these low-order degeneracies. 

As an example we consider a model for the expansion in which the free parameters are 
Qm, and the parameters wq and Wa for the dark energy equation of state w = wq + Wa{l — a). Note 
that Dde = 1 — — ^k- We will also quote the uncertainty in Wp, the equation of state at the 

“pivot redshift” Zp = l/op — 1 for which the model w = Wp + (op — a)wa has uncorrelated Fisher 
errors in Wp and Wa- We consider the constraints that arise from our canonical survey of weak 
lensing cross-correlations and photo-z BAO data. We presume a hducial proxy-mass correlation 
coefficient of r = 0.7. Table 1 lists the Fisher-matrix estimates of the standard deviations on these 
four parameters. Figure 3 plots error ellipses in the WQ—Wa plane after marginalization over the 
other two variables. 

We find first that the WL CCC alone retains a strong degeneracy between and the other 
parameters, with large errors of 0.36, 0.60, and 2.0 on 17^, and Wa, respectively, but with 
correlation coefficients > 0.99 among all of these. [These are for full-sky surveys; we will omit the 
factor /gky"^ brevity here.] The degeneracy is broken either by assuming a flat Universe, by using 
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Fig. 3.— Error ellipses in the WQ-Wa plane for parametric models of dark energy that have been marginalized 
over flm and fife. These are 68% CL regions (Ay^ = 2.3), for full-sky coverage; errors scale as /sky^^- The 
large ellipse is for WL CCC information only. The middle shaded ellipse (red) adds transverse baryon 
acoustic oscillation data from the photo-z information of the WL survey. The innermost (green) shaded 
region includes a measurement of the distance to recombination from expected CMB data. The two dashed 
ellipses allow for systematic shear calibration errors of 3% per redshift bin, which is in the “self-calibration” 
regime. The outer ellipse is CCC-I-BAO, the inner adds the recombination distance. 


the CMB distance to recombination, or by incorporating the transverse BAO information over the 
same sky fraction. In the last case, has uncertainty 0.007, which as expected is a much tighter 
constraint than the model-independent one above. We also note that the Wp, Wa uncertainties of 
(0.013, 0.057) are only slightly reduced by assuming flatness, although their Zp shifts. 

Inclusion of an accurate (0.1%) angular-diameter distance to recombination from the CMB, as 
might be expected from Planck, drops uncertainties on by a factor 3, and a factor 2 for Wp. 
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4.1. Scaling with Survey Characteristics 

The BJ04 analysis asserts that CCC constraints improve without bound as ZHS did 

not concur, and we would now agree. Since Stebbins (private communication) proves that G must 
be non-singular in the continuum limit (and we have constructed its inverse), then the matrix 
K = N -|- GCuG^ will remain non-singular as the shape noise N ^ 0, meaning that the Fisher 
information in Equation (45) remains finite. As a practical matter, we find that in the range of n 
accessible to planned observations, CCC Fisher information scales only slightly more slowly than 
^- 1 / 2 ^ There is no simple scaling once we combine the CCC information with that from BAO. The 
constraints are nearly independent of the choice of 6^1, again because the shape noise dominates 
sample variance. 

The BJ04 derivation also leads to CCC-only parametric uncertainties that scale nearly inversely 
with the fiducial correlation coefficient r. In the present analysis we see that such scaling arises if 
we consider only the effect of r upon the Cg factor of Equation (46), and if we drop the left-hand 
term entirely. But if the effect of r upon the Cu terms in the noise matrix K are also considered, 
then the dependence is more complex. We find that the CCC-only parametric constraints on flat 
Universes scale more weakly than inversely with r over the range 0.3-1.0. The dependence is yet 
weaker once we include the BAO information. 


4.2. Effects of Calibration Errors 

Eigure 4 plots the degradation of the parametric cosmology constraints as we weaken the prior 
on the calibration factors /*. There are three regimes: when the calibration is known a priori to 
better than there is no significant degradation due to calibration errors. Over the 

range < cj/ < 0.01/^]^^'^^, the constraints continually degrade. Calibration uncertainties 

above 0-01/sky^ > however, cause no further degradation in cosmological constraints, as the CCC 
data solve for the systematic factors in a “self-calibration” regime. The self-calibration regime has 
parameter errors that are 1.6-2 times larger than the no-systematic regime. We conclude that 
lowering the calibration errors from 0.01 to 0.001 is equivalent to increasing the survey area by 
a factor of 2.5-4, but interesting cosmological constraints are possible even with poor a priori 
calibration systematics. 

We find a similar scaling of uncertainties with calibration systematics for a shallower, ground- 
based survey. 


4.3. Comparison to Previous Work 

We alter the parameters of our fiducial survey to match those of ZHS, namely we set /gky = 0.1, 
zq = 1.0, = 0.3/\/2, Hfc = 0, and we use an equation of state w = wq + w'z instead of the Wa 
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Fig. 4.— The l-a uncertainties in each of the labelled cosmological parameters in a CCC+BAO survey are 
plotted vs the accuracy of our prior knowledge of the shear calibration factor fi in each redshift bin. We see 

that calibration errors smaller than « are harmless, and there is a “self-calibration” regime for 

_ /2 

CT/ > O.Ol/gjjy , in which the CCC data solve for the calibration factors to an accuracy that is better than 
the prior. The degradation from self-calibration is mild. Both axes may be scaled jointly with ■ 


form. The ZHS analysis also allows for an arbitrary distribution of mass redshift within a shell, and 
accounts for the signal attenuation due to photo-z errors. We compare to their = 0.01 case. For 
a fixed-rc model with only CCC information, we obtain standard deviations on {Qrn,u!) of (0.028, 
0.046) for r = 0.7, roughly 2x smaller than ZHS derive. We must lower r to 0.25 to obtain similar 
results. They also derive errors on (wq, w') when marginalized over with a prior of 0.03. In this 
case our r = 0.7 constraints of (0.03,0.11) compare to their (0.07,0.09)—similar in size but different 
in shape. 

The agreement to within factor two must be considered adequate at this time. While there 
remain substantial differences in the analytical approach to the CCC method, the discrepancies 
could easily be attributed to different assumptions about the correlation between the mass and 
the proxy field. ZHS, as well as Hu & Jain (2004), take the galaxy distribution itself as the proxy 




- 21 - 


field, and assume that the galaxies populate the mass via a pure Poisson process. In this case the 
correlation coefficient between multipole I of the proxy and mass fields at distance Da is 


f 

IdN , , J 


yp(‘/DA\ 



(54) 


where P{k) is the 3-dimensional mass power spectrum. Thus when nP < 1, r drops below the 
0.7 value of our fiducial model. This occurs at I > 10^ in our fiducial model for 0.5 < z < 2, 
whereas our fiducial model makes use of shear power to higher 1. Were we to adopt a Poisson 
model our parametric dark energy constraints would weaken. It is not clear, however, that the 
Poisson model is appropriate once we approach the length scales of group- and galaxy-scale halos, 
which are currently suspected of harboring a single central galaxy with near-unity probability, plus 
satellites. The galaxy-mass covariance on small scales may thus be a good deal stronger than 
Poisson models predict. This is surely a subject for further research. 

We may also compare to the CCC-only results with the BJ04 method by reverting to the 
parameters assumed therein for the “SNAP” case; /sky = 0.025, r = 0.8, = 0.15, 12^ = 0, 

and a prior of 0.03 on Qm- We find a pivot redshift Zp and uncertainty cr{wp) that are the same 
in both cases, but Wa is 35% lower with the present Fisher matrix. This could be a result of 
some of the simplifications made by BJ04 in construction of the covariance matrix for the Xgi in 
Equation (30). Alternatively it could derive from our present assumption that the convergence, 
mass, and proxy field values have a multivariate Gaussian distribution, while BJ04 have the less 
restrictive assumption that the Xgi have a Gaussian distribution. 

A close comparison with Hu & Jain (2004) is well beyond our present scope, because these 
authors include a parametric model for the bias function and limit the cross-correlation information 
to multipoles I < 3000. In principle both differences can be accommodated in the present formalism, 
by switching to a spherical-harmonic decomposition, and by projecting parametric bias constraints 
onto the Fisher matrix before marginalizing over the Bi. A WL GGG survey reduces 

the errors on InH to < apart from a quadratic degeneracy with distance. Hence the a 

priori theoretical constraints would need to be stronger than this in order to improve dark energy 
measurements. We also note that Hu &: Jain take a Poisson model for the mass-galaxy correlation, 
and assume a density of lens galaxy samples that is 2 orders of magnitude lower than assumed 
here, which we believe accounts for the cross-correlation constraints being much weaker in that 
work than in the present calculations. 

A short summary of this subsection is that there are currently multiple formalisms for analysis 
of CCC data, and they currently differ by up to a factor 2 in Fisher uncertainties on parametric 
dark-energy models. Further work is certainly required in order to understand the nature of the 
CCC constraints, and to understand our ability to reconstruct the mass distribution from the galaxy 
data. 
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5. Discussion 

Reparameterizing the observable galaxies by their (dimensionless) angular diameter distance 
d rather than redshift z makes it clear how the WL CCC method can provide model-independent 
geometric constraints on This methodology also reveals that CCC data is degenerate under 
alterations to In 6 and Ind by quadratic functions of d, if is free. Apart from three exact 
degeneracies, the cross-correlation strength, being a joint function of the lens and source distances, 
is very efficient at producing decoupled estimates of 0^, the distances di and the bias factors Bi, 
and can even determine shear calibration factors fi for each source plane with modest degradations 
of factor < 2 in cosmological accuracy. 

With sufficiently accurate photometric redshifts, the same survey data used for weak lensing 
may be used to determine the transverse BAO scale. This gives additional model-independent con¬ 
straints on d that are of lower precision than the CCC data, but are completely free of degeneracy. 
Such a combined CCC-BAO survey should yield uncertainties of ss 0.04/^^,^ on We reiterate 
that such constraints are completely independent of any assumptions about the matter-energy con¬ 
tent of the Universe, any biases in photometric redshifts, or in fact any alterations to the Friedmann 
equations or the deflection equations for light. They merely require that the Robertson-Walker met¬ 
ric be applicable to our Universe. Given the lack of viable dark-energy theories, it seems prudent 
to seek cosmological information that remains valid even if the acceleration is attributable to an 
alteration of General Relativity rather than a previously unnoticed stress-energy contribution. 

In the spirit of conservatism, we should note that the RW metric may not be sufficiently 
accurate. The dumpiness of the matter-energy distribution may invalidate our adoption of the 
filled-beam angular diameter distance, and a more sophisticated treatment may be required (Holz 
& Linder 2004). It is not clear to what extent small-scale dumpiness can influence the apparent 
shear of galaxy-scale images as surveyed across the entire sky—this is quite a different regime from 
studies of strongly lensed quasars or supernovae for which the Dyer-Roeder distances have been 
most carefully studied. 

The curvature measurement is limited by the ability to break the WL CCC degeneracy, i.e. by 
the accuracy and redshift span of the BAO or SN measurements of angular-diameter distances. Our 
baseline constraint assumes photo-z BAO information for 0 < z < 3 over the full sky. A massive 
spectroscopic BAO survey would offer substantial improvement, but Type la supernova studies 
would not offer substantial improvement over full-sky photo-z BAO unless systematic uncertainties 
in SNIa peak magnitudes could be brought well below 1% and the redshift range extended beyond 
z = 2. Inclusion of lensing and acoustic-scale information from the epochs of recombination or 
reionization would help significantly. 

The parametric dark-energy constraints derived from the CCC-I-BAO information are more 
fluid, as our derivation has made somewhat arbitrary assumptions about the correlation between 
the true mass and the proxy mass derived from the galaxy distribution. The simplifying assumptions 
about Gaussianity, uncorrelated mass shells, and uncorrelated lines of sight should also be improved 
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upon. The spherical-harmonic-based Fisher analyses of Hu & Jain (2004), Song & Knox (2004), 
and ZHS are in this respect superior to our treatment of uncorrelated lines of sight. We note 
however that the Fisher matrix from Equation (45) is applicable to multipole coefficients. More 
work is also required to examine the impact of finite photo-z precision on the CCC constraints, as 
per ZHS. 

A weak leasing survey would produce other kinds of information which we have not used in 
this analysis. Power-spectrum (Hu 1999; Refregier et al. 2004) and bispectrum (Takada & Jain 

2004) tomography, and WL cluster counts (Weinberg & Kamionkowski 2003; Marian, & Bernstein 

2005) constrain dark energy once a model for the growth of structure G{z) is specified. One could 
augment our current 0^ — di Fisher matrix with growth factors G{di) to produce a non-parametric 
information estimate, which could then be projected onto the parameters of specific models, even 
those which predict changes to the usual equation for linear growth of perturbations in an expanding 
Universe. 

The method proposed herein for metric constraints on curvature, and for strong parametric 
constraints on cosmology, will require substantial practical advances before implementation, plus 
additional theoretical work. Progress in high-accuracy shape and photo-z measurement is needed, 
as is research into intrinsic galaxy-shape correlations (Hirata & Seljak 2004), non-linearities in 
lensing (White 2005), and methods of reconstructing the matter distribution from a galaxy distri¬ 
bution. Nonetheless the WL CCC method shares with the BAO method the ability to generate 
extremely robust and model-independent cosmological constraints, particularly when they are used 
in combination. 
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Table 1. CCC+BAO Constraints on Cosmological Parameters 


Photo-z Error 

Calibration Prior*^ 

a{Qk) 

Parameter Errors'^ 
cr(Dm) cr{wp) 

cr{wa) 

CCC only 


0.36 

0.60 

0.018 

2.0 

CCC only 



0.015 

0.017 

0.081 

0.04 


0.0066 

0.0045 

0.013 

0.057 

0.04 



0.0038 

0.010 

0.052 

0.04 

0.03 

0.013 

0.007 

0.044 

0.100 

0.04 

0.03 


0.007 

0.026 

0.083 


®'A11 uncertainties are for the fiducial survey described in the text, with 
full-sky coverage, but scale with as long as the calibration prior is 

similarly scaled. Dashes indicate quantities that are held fixed. Each error 
assumes marginalization over all the other parameters. 








